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FINAL TECHNICAL REPORT 


CONTAINERLESS PROCESSING TECHNOLOGY ANALYSIS 


NASA CONTRACT NAS8-33731 




by J. E. Rush, Principal Investigator 


I. INTRODUCTION 

This report covers research on acoustic levitation, air-jet 

j* J 

levitation, ■ and heat transfer from molten samples. Although the work on 
these topics was not completely sequential, they are separated in this 
report for clarity. The performance period was May 8, 1980 to August 
31, 1982. 

The thrust of this research was toward obtaining a better 
understanding and improving the quality of containerless processing 
systems of interest to Marshall Space Flight Center (MSFC) . Such 
systems have application to the study and processing of materials in 
situations in which contact with a container must be avoided, and have 
potential application in both ground-based and orbiting laboratories. 
An overview of the general subject of containerless processing is given 
by Naumann and Herring (1). Typical applications of the systems studied 
here are in the development and study of glasses from materials which 
normally crystallize upon cooling. 

In addition to the reports which have resulted from this work and 
which are identified in the following sections, the PI gave 
presentations on this work at the Space Sciences Laboratory, MSFC; a 



2 


Materials Research. Society Symposium (see Appendix B) ; the University of 
the South Sigma Xi Chapter; the University of Alabama (Tuscaloosa); and 
The University of Alabama in Huntsville. 

II. ACOUSTIC LEVITATION 

A. Background 

A single-axis acoustic levitator has been under study and 
development for several years under NASA sponsorship, the primary 
contractor being Intersonics, Inc. The work of this report has been 
primarily in support of, and in addition to, that of Intersonics. The 
Principal Investigator (PI) is indebted to Dr. Roy TThymark and Dr. 
Charles Rey of Intersonics for their cooperation in sharing information 
during the course of this research. 

The acoustic levitator utilizes a high- intensity standing-wave 
sound field to levitate or position small (1 cm or less) objects. The 
basic device was designed by St. Clair (2) in 1940 for a completely 
different purpose. Current levitator designs are available In reports 
from Intersonics, Inc. The basic theory for employing high- intensity 
sound sources for levitation was worked out by King (3) in 1934, based 
on earlier work by Rayleigh. 

The PI had done work related' to the levitator prior to the award of 
this contract, first as a NASA/ASEE Sunimer Faculty Fellow at MSEC (1979) 
and later as a physicist on a related contract. The fellowship involved 
experimental work relative to the causes of problems in spot-heating 
samples in the levitator at 1 g, and is summarized in a final report (4) 
and in the proceedings of a symposium (5). Earlier work on the 
levitator was reported by Oran, et al. (6), and by tThymark (7). 
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The work following the fellowship- was directed at obtaining a 
better theoretical understanding of the failure of spot heating. It 
turned out that quite a bit of related research had been done in other 
contexts, primarily involving the effect of acoustic fields on heat 
transfer from cylinders. Most of the pertinent references are contained 
in a review article by Richardson (8)'. The significant effect noted by 
these authors is a dramatic increase in heat transfer from cylinders in 
air when placed in an acoustic field of 140 db or greater. 

I'/hile several explanations for the enhanced heat transfer rate have 
been proposed and the definitive answer has not been given, it is clear 
(9) that the acoustic field causes the laminar flow of gravitational 
convection to become turbulent, and this turbulence is correlated with 
the drastically increased rate of heat, transfer. It was this turbulence 
which caused spot-heated samples to be lost from the 170 db acoustic 
field in the laboratory. The turbulence, however, is not manifest in 
acoustic fields with' unheated samples. The fact that neither heating 
nor high-intensity sound, independently, produce significant turbulence, 
while the combination does, is not particularly surprising because of 
the nonlinearity of the equations at the intensities employed for 
levitation. Indeed, the nonlinearity is basic to the levitation process 
(see the equation for the mean sound-pressure level derived by King, 
Ref. 3, p. 215). 

A summary of the published results on heat transfer in acoustic 
fields, including a brief analysis of each paper., is available from the 
PI. Also available are 16 ram motion-picture films showing the behavior 
of samples at 1 g, including the effect of heating-, and plots of dc 
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(time averaged) pressure, ac pressure, and temperature in the sound 
field, with and xjithout spot-heated samples. 

B . Further Studies of Heat Transfer in Acoustic Fields 

Because no definitive conclusions about the mechanisms for heat 
transfer in acoustic fields had been reached in the published 
literature, and because the parameter space involved in our levitation 
work had not been explored by other authors, it was considered advisable 
to carry out further studies of the phenomenon. The apparatus consisted 
of a very small cylindrical solenoid of Kanthal wire coated with a 
nonconducting material to form a l.-O cm long cylinder with a diameter of 
0.25 cm. The cylinder (several similar ones were used) was attached to 
nichrome wire and rigidly supported in the sound field at a velocity 
antinode in the standing waves. 

Detailed measurements were made of the heat transfer rate in 
thermal equilibrium for cylinder temperatures up to 700°C and sound 
pressure levels- (SPL) up to 170 dh. Higher SPL values were available 
from the apparatus but could not be accurately measured. An estimated 
SPL of 178 db was observed. 

Further details of these measurements, results, and comparison with 
published papers, are contained in an interim report which is made 
Appendix A of this final report. Note that Appendix A has a 
self-contained list of references. 

C . Analysis of Behavior of Sample in Acoustic Field at Low G 

At the request of R. Naumann of MSFC, the PI assisted in the 
analysis of results of a single-axis levitator experiment on the 
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SPAR-VIII flight (^ace ^Processing _^plications Rocket). The experiment 
was No. 74-42/2, flown on November 18, 1980. A preliminary analysis had 
already been carried out by C. Rey of Intersonics, Inc. The additional 
analysis was done by the PI and by C. F. Schafer of MSFC, with the 
assistance of R. L. Holland of MSFC. The results are presented in a 
NASA Technical Memorandum (10). 

In addition to the experiment analysis and recommendations for 
improvement in future SPAR flights, the memorandum contains unpublished 
results on the radial properties of the sound field in an acoustic 
levitator, an extension of the work of King (3) , done by the PI. 
Experimental vrork on acoustic streaming in the sound field of the 
levitator was also, carried out by the PI in' collaboration with C. F. 
Schafer, and with W. K. Stephens of UAH. This work was not completed 
because of more pressing matters, but preliminary results, including 
photographs and 16 mm film, are available from the PI. 

Copies of the NASA Technical Memorandum are available from the PI, 
as well as from the^ customary sources. 


D. Other Work on Acoustic Levitation 

In connection with the work discussed in Sections A-C, several 
related measurements were made, some of which are described in Appendix 

A. 

A frequency analysis of the acoustic driver and the standing-wave 
field was made to facilitate an understanding of the nature of the 
field, and the results are- in Appendix A. A study was also made of the 
lateral force on a small sphere, in connection with the analysis of 
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SPAR-VIII. The results of this study are contained in laboratory 
notebooks available from the PI. 

III. AIR- JET LEVITATION 

A. Background 

The technical monitor for the initial work on acoustic levitation 
was W. A. Oran of MSEC. Following Dr. Oran's assignment to NASA 
Headquarters in September 1980, the technical monitor has been Dr. E. C. 
Ethridge of MSFC, All work by the PI on the air-jet levitator was done 
under Dr. Ethridge. The work was initially begun at the request of Dr. 
R. J. Naumahn of MSFC. 

A constricted-tube gas-flow levitator was developed at MSFC by 
Berge, Oran, and Theiss (11). The device consists of a quartz tube with 
a constriction introduced by melting and stretching (tubes were prepared 
by R. Smith, glassblower at UAH); a source of compressed air or gas; and 
a furnace for heating the air before it passes through the tube. By 
means of this device, Berge, et al. , were able to levitate a spherical 
sample with the tube upright or inverted, and to heat the sample to 
temperatures greater than 1300° C, 

The only results in the open literature directly related to this 
device were by Schmidt and Springer (12), involving a sphere in a 
diffuser, and they are not particularly useful for the air-jet (or 
gas-flow) levitator. There is, however, an extensive literature on 
spheres in tubes (see references in Appendix B) . 

The PI was asked to predict the behavior of samples in the 
levitator in a low-g environment. 
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B. Behavior of Samples in Air-Jet Levitator in Low G 

The experimental parameters involved in understanding the behavior 
of the air-jet levitator are: 

(a) flow rate of air or gas 

(b) pressures in tube 

(c) relative diameters of tube and constriction of sample 

(d) angle and shape of constriction 

(e) net force on sample at different positions in tube 

(f) viscosity of air or gas 

In order to' predict the behavior' of samples at low g, a study was made 
in which (b) and (e) were measured as (a), (c), and (d) were varied. 
The effect of (f) can be calculated from the Reynolds number. 

The results of these measurements were presented at the Materials 
Research Society Symposium, "Materials Processing in the Reduced Gravity 
Environment of Space," Boston, November 16-18, 1981, with W. K. Stephens 
of UAH and E. C. Ethridge of MSEC. They are given in Appendix B, which 
is a copy of the paper which was presented. The conclusion is that the 
levitator should work as a positioning device at low g, but that care 
must be taken to regulate flow rates so that the molten sample never 
touches the tube. The development of the device for low g will thus 
require a careful study of flow rates vs. temperature, in a low gravity 
environment, .to maintain sample stability. 

C. Behavior of Samples in the Laboratory 

S.tudies’ were also made of sample behavior in the laboratory at 1 g, 
and some of the results are contained in Appendix B. It is possible to 
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do a theoretical analysis of the air-jet levitator for laminar flow with 
a precisely defined geometry, and such an analysis would be quite 
similar to flow-meter theory. However, the flow is not laminar in, the 
actual levitation experiments, and the academic interest of such an 
analysis hardly seems appropriate in the context of developing working 
devices . 

The furnace arrangement built by Berge, Oran, and Theiss was 
redesigned and rebuilt by the PI with the help of R. Eakes of UAH, and 
is available for continued work at 1 g, 

IV. HEAT TRANSFER FROM A MOLTEN SAMPLE 

A. Background 

The basic problem involved in this part of the project was to 
develop a computer model for the temperature distribution, as a function 
of position and time, in a small spherical drop of material which is 
cooling by conduction (in the material) and radiation. The material is 
neither opaque nor fully transparent. This problem is related to the 
processing of samples in drop tubes, on stings, etc. 

The problem described above has never been solved. However, there 
are several simpler, related problems that have been solved approxi- 
mately. The subject is discussed generally in a review article by 
Viskanta and' Anderson (13). The basic problems of radiative transfer 
are presented by Chandrasekhar (14). 

Some steady-state problems for which approximate solutions are 
available are 

1) a gray spherical shell with no conduction (15-17); 

2) a special nongray spherical shell with no conduction (18) ; 
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3) a gray spherical shell with radiation and conduction (19, 20); 

4) conduction with no radiation (21) . 

Time-dependent problems which have been considered are 

5) cooling of a spherical gray gas with no conduction (22, 23); 

6) conduction in a sphere with no radiation (24) . 

B. Theoretical Model 

In order to simplify use of the word processor for this report, we 
modify the conventional notation somewhat. Some pertinent variables in 
the problem are 

u - energy density 

^ - heat flux vector 

t - time 

D — density 

c - specific heat 

T - temperature (absolute) 

^ “ frequency of radiation 
k - thermal conductivity 

K - linear radiative absorptivity (depends on f) 
n - index of refraction (depends on f) 

The energy conservation equation is then 

du/dt - “ div q 

where, for this section, d implies partial differentiation. Then from 


q = - k grad T + F 
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and u = DcT 
we get 

Dc dT/dt = div (k grad T) - div 
which, with spherical symmetry, becomes 

Dc dT/dt = r“^ ^ (kr^dT/dr - r^F) 

where F is the radial component of and r is the radial coordinate. 

The problem of the distribution of radiant energy for a given 
temperature distribution was solved by Chandrasekhar (14). The term div 
F^ can be represented as 

y * 00 9 , 

K(f)[4Trn (f)l(f) - G(f)]df 
o 

where 1(f) is the blackbody distribution function and G(f) is the 
contribution from absorption in each small volume of material. The 
expression for div ^ in terms of temperature is given by Viskanta and 
Lall (22) in terms of exponential integral functions. 

For the simplest case, we assume that c, k, K, and n are constant. 
Then we can change to dimensionless variables 
X = r/r 

o 

w = xT/T 

o 

t* = t/(Dc/KBT ^) 
o 

N = kK/BT ^ 
o 

s = t* N/(Kr )^, 

0 

where is the spherical boundary, T^ is the initial temperature, and B 
is the Stefan-Boltzmann constant. Then we get 
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du/ds = d^u/dx^ - H 

where H is an integral representing the flux, which depends on position 
and temperature. 

If we define 

y = Kr 


so that X = y/y 

o 

we find 

H = y(4E'-G') 
where E' = (T/T 

o 

and G' = G/BT ^ 

o 

with- G the integral over frequency of G(f). Viskanta and Lall get for 
G, 


G = (2/y) J^^E' i\y~y' \) 
o 

- \(y + y')3 y'dy' 

where E^^ is the exponential integral function. 

The basic equation is thus an integro-diff erential equation of the 
parabolic type. We have explored a variety of possibilities for 
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approximate solutions to this equation, and finally settled on a direct 
computer approach using finite differences and numerical integration. 

The most significant approximations which have been made are taking 
k independent of T and taking K independent of f. Relaxing the first 
restriction would make the calculation slightly more complicated, but 
the equation is already strongly nonlinear because of the flux integral, 
which is not small, so it would be a minor complication. Relaxing the 
second restriction is quite a serious matter. 

IThether or not the second approximation is reasonable depends on 
the shape of the absorption curve for a given material in the vicinity 
of the blackbody peak. The position of the peak follows from Wien's 
displacement law 

SL T = 2898 pm-K 
max 

for wavelength St,. For the temperature range 2000 K to 3000 K, this 

gives Si from 1.5 to 1 microns. We have as yet no data on K for 
max 

liquid alumina, but for solid alumina the transmittance is essentially 
uniform from 0.2 to 5 microns. To this extent the approximation of 
constant K for alumina ^ good one. For each additional 

substance of interest, the approximation must be reconsidered. 

From an analysis of energy transfer at the spherical surface, it 
would appear that the appropriate boundary condition is 

dT/dr = 0 at r . 


However, thiS' gives 
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dw/dx = T/T at r , 

0 o 

which is impracticable. A relatively minor modification is 
dw/dx = 0 at r 

o 

which gives 

dT/dr = - T/r at r . 

o o 

Starting in thermal equilibrium at t = 0, we also expect 

T = T for r between 0 and r , but this gives a null solution, 
o o 

so we have chosen 

T = exp (-aT) 
for a near zero. We also have 

w = 0 for x = 0, all t. 

The difference equation is given, e.g., by Ames (24). The computer 
program (in Fortran V) is given in Appendix C. We have incorporated the 
necessary condition for convergence (see Ref. 24, p. 323). 

The general result of the computer calculation is what one would 
expect. The temperature at any point decreases with time. At a given 
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time, the T vs. r curve looks roughly like a decreasing exponential 
function. 

For sample numerical results, we used the following values: 

melting point of ~ 2323 K 

T = 2500 K 

o 

3 

density D = 4.8 g/cm (Ref. 25) 
heat capacity c = 197 J/mole • K (Ref. 25) 

thermal conductivity k = 0.065 W/cm • K (Ref. 26) 

optical absorptivity K = 0.2/cm (Ref. 27) 

We find 

N = 0.13 
t =1.5xl0^t 

We divided x into 50 equal segments. To satisfy the convergence 

-4 

requirement, this means that a good value for an Increment of s is 10 
or less. Then we get 

t = (140 sec/cm^) r ^ 

o 

For r^ = 1 mm, incrementing t 50 times, we get a total time of 72 sec, 
which should be quite adequate. 

Since the computer program involves dimensionless variables, a few 
runs for typical values of y^ are sufficient to cover all cases of 
interest with different substances. The results have been made 
available to the technical monitor at MSFC. 
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The PI is continuing to search for improved experimental data on 
substances of interest to MSFC and to study ways of expanding the 
applicability of the model. Any additional results will be made 
available to appropriate personnel at MSFC. It is intended that the 
results of this xrork be published. 
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Interim Report on Contract NAS 8-33731 
J , E . Rush 

This Is a much-expanded version of the first quarterly report on 
subject contract, covering work for the period June 2 - August 7 (no 
work done on the contract for May 8 - June 1) . All the work involved 
the two acoustic-levitation systems at SSL/MSFC and included: 

(a) instrument analysis 

(b) experimental study of heat transfer from cylinders and spheres, 
and 

(c) calculations of relevant physical quantities and comparison 
with published data. 

Detailed data are contained in my laboratory notebook No. 2. 

A. Instrument Analysis 

Studies were done of the stability of the Intersonics (Model 15-lA) 
feedback amplifier and driver system (St. Clair generator). It was 
determined that a cooling fan directed at the amplifier could reduce 
the downward power drift of - 4 db in 8 min to ~0.5 db, and reduce 
the rapid fluctuations from + 0.3 db to + 

During the course of these studies, an instability in the driver 
reflector system was noted. This consists of a rapid fluctuation of 
~ 4 db when the system is tuned to resonance, at ~178 db. It occurs 
with a concave driver and large flat reflector or a flat driver and 
concave reflector. It was eliminated by use of a flat driver and flat 
reflector or a concave driver and small flat reflector. Studies with 
oscilloscope and wave analyzer showed that the instability involves 
alternate coupling and decoupling to the first and second harmonics of 
the fundamental vibrating frequency (~ 15 k Hz) . 
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Because the B & K microphones were registering 178 to 180 db for 
the system in resonance, and they are inaccurate beyond 170 db, an effort 
was made to obtain an alternate source for measurement of high sound 
pressure levels (SPL) . With the help of John Theiss of MSFC, I obtained 
a strain-gauge pressure transducer, a sensor amplifier module (SAM-1) 
to accompany it, an oscillograph to record the signals, and a deadweight 
tester and 1 kHz standard source to calibrate it. The transducer 
was rated from 0 to 5 psi (184 db) and 0 to 40 kHz . 

tflien the dc signal was calibrated with the deadweight tester, 
the transducer agreed well with the microphones at 1 kHz and 160 db, 
using the standard source. Unfortunately, it gave readings about 4 db 
higher than the microphones when placed in the field of the acoustic ■ 
system, at SPL values from 155 to 180 db. The signal from the trans- 
ducer was also quite different from the microphones as seen on the 
oscilloscope and the wave analyzer. Since the B & K microphones are 
knoTO to be accurate below 170 db, I was forced to discard the strain- 
gauge transducer as a reliable measuring instrument for the acoustic 
system, and be left with no accurate measuring device for SPL values 
greater than 170 db. Because of instabilities at higher SPL values, 
this limitation may not be a practical problem. 

During the above calibration process, I did wave analysis studies 
of the drivers with no reflector and clearly observed the sawtooth effect 
predicted for high SPL values . 

B. Experimental Study of Heat Transfer 

As a continuation of studies begun in the summer of 1979 and 
continued during the period September 1979 to May 1980, I constructed 
several heating elements and measured their heat transfer characteristics 
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In the acoustic field. The configurations were as given in Table 1. 


TABLE 1 

Coil Configurations 


Shape 

No . Turns 

Length 

Width 

sphere (approx.) 

5 

0.4 cm 

0.4 cm 

cylinder 

20 

1.0 cm 

0.3 cm 

cylinder 

20 

1.0 cm 

0.25 cm 


In each case the coil was wound with B & S #28 Nichrome wire 
and spot-welded to leads of B & S # 20 Kanthal wire. The coil was 
dipped in Sauereisen resistor cement, shaped, and allowed to dry with 
an iron-constant an thermocouple, made of 0.004 in. wire, imbedded 
approximately in the center. The leads were then inserted into a rigid 
ceramic holder which was taped to a rod which was mounted on a ring 
stand. Power was provided by a very steady (+ 0.003 A) d.c. power supply. 
The thermocouple leads were connected to an Omega 2168A digital thermom- 
eter which showed integer degrees. 

The initial heat -transfer measurements were made in a sound field 
with SPL recorded at ~ 178 db by B & K 1/8" microphone (see notes 
about SPL measurement in Part A) , using the Intersonics Model 15-lA 
systems. A large flat aluminum reflector tras placed about 3.0 cm from 
the edge of a concave aluminum driver. The reflector was then raised 
slightly ( ~ 2 mm) so that the system was just far enough from resonance 
to avoid the large pressure fluctuations described in Part A. Several 
sets of data were taken with the sphere and the second cylinder, the 
thermocouple in the first cylinder having failed. Measurements were 
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also made of heat transfer from the sphere with no sound field, with 
no pieces of apparatus near enough to influence convective heat trans- 
fer. In the sound field, the coil was placed at the location of the 
first minimum below the highest maximum below the reflector, putting 
it approximately halfway between driver and reflector. This location 
should give maximum velocity of the a.c. field and maximum cooling. 

(This fact had been checked in earlier measurements . See laboratory 
notebook No. 1.) 

In order to determine the effect of radiation and conduction of 
heat through the leads, the emissivity of the sphere was measured with 
an IRCON radiometer calibrated to a standard blackbody source; and 
thermocouples were attached at each end of a Kanthal lead with Sauer- 
eisen to determine the temperature gradient for various temperatures 
of the sphere. (These connections and measurements were independent 
of the convection measurements . ) 

The current to the coils was measured with a Hewlett Packard digi- 
tal ammeter for all values _< 2.00 A, and the voltage across the coil 
was estimated by sliding probes connected to a digital voltmeter along 
the Kanthal leads and extrapolating to the limiting value at the points 
of attachment of the coil. Current values greater than 2.00 A were 
gotten bysetting the ammeter on the power supply toa corrected value estimated from 
lower current readings withthe digital anmeter. (Only a zero-pointcorrectionwasneeded) . 

Following the measurements at - 178 db, measurements were made at 
170, 160, 150, 140, and 130 db. Those at 170 db were made with the 
system configuration noted above, with the reflector adjusted for lower 
SPL. The remainder were made on the second amplifier-driven system, 
using an ALTEC 9440 A, 800 watt amplifier with oscillator and magnet 
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circuits constructed at MSFC. The second system was used to allow 
stability at lower power. To arrive at 130 db, the driving cylinder 
had to be firmly clamped. At 130, 140, and 150 db, the signal read 
by the microphone was an almost pure sine wave at 15.2 to 15.5 k Hz 
(depending on driving cylinder) . At 160 db it was slightly distended 
toward a sawtooth. At 170 and 178 db it was clearly distorted. All 
SPL readings were taken with a 1/8" B & K microphone connected to 
a B & K Type 2607 measuring amplifier. The microphone was attached 
to a cathetometer so that it could be easily moved in and out of the 
sound field. All SPL readings were made at the maximum SPL point, and 
the microphone was removed before temperature readings were taken. 

All readings were rms pressure . The shape of the signal was 

determined by feeding the output of the B & K measuring amplifier 
into an oscilloscope. 

Comparison with a good sine wave was easily made by inserting a 
22.5 k Hz low-pass filter (included in the B & K unit) or by passing 
the signal through a wave analyzer and looking at the output of the 
wave analyzer on the oscilloscope. A dual trace unit on the oscillo- 
scope was used for easy comparison. 

C. Heat-Transfer Calculations 

Most of the calculations were done using an HP 9835 A desktop 
computer located in SSL Room 215, by means of five programs - labeled 
SPL, QTRANS, QCOR, QMEAS, and QCAL, which I wrote and which are stored 
on an HP cassette. The programming language is BASIC and the pro- 
grams can be easily read and interepted by anyone familiar with FORTRAN. 
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1. Free Convection 


For free convection, Yuge (1) has obtained an empirical 
formula for spheres, 

Nu = 2 + 0.392 1 < Gr < 10^ 

where Nu is the Nusselt number and Gr is the Grashof number. 
His defining equations are 



m 


and Gr = g D — 2 1 ■ 

m 

Here D is the diameter, h the heat-transfer coefficient, k 
the thermal conductivity of air, g the acceleration due to grav- 
ity, A t the difference between surface and ambient temperature, 

V the kinematic viscosity of air, and T the absolute ambient 

a 

temperature. The subscript m indicates the value at the mean 
between surface and ambient temperatures (the film temperature) 
and h is defined by 

^ AAt 

where Q is the heat-transfer rate and A the surface area 
(= IT D^) . 

A complete list of symbols with definitions is given at the 
end of this report. All Grashof numbers for my data were between 
20 and 200, thus falling within Yuge's limits. 

The free-convection measurements which I made agree with 
Yuge’s formula to within 16% at all temperatures up to 550 "C, 
even though his data only went to yO^C. In the comparison, the 
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radiation and conduction losses were subtracted from the measured 
power input to the coil. The conduction loss was obtained from 


Q 


c 


2 k 


u) 


A 

u) 


Atoj 

L 


CO 

where A is the cross-sectional area of the Kanthal wire, K 

ul U3 

the thermal conductivity, and At^ the temperature difference 

over the length L^. This does not allow for convective cooling 

of lead wires, and so overestimates the convection loss of the 

sphere. This fact coincides with the data, which show convection 

losses progressively larger than those calculated from Yuge’s 

formula as the temperature is increased. At 550° C, the calculated 


radiation and conduction losses accounted for 19% and 18%, respect- 
ively, of the measured power input to the coil. 


2. Convection in the Sound Field 

labile there do not appear to be any published data on heat 
transfer from spheres in an acoustic field, there are several 
published papers dealing with heat transfer from cylinders, either 
in a sound field (2-4, 6-12) or vibrated (5,7). For sufficiently 
long wavelengths, there should be no distinction between the 
effects of sound and vibration (5) . Ford and Peebles (7) determi- 
ned that the effects were indistinguishable for X/D ^ 12, where 
X is the effective wavelength and D is the cylinder diameter. 
The X/D value for my data is 9, so some distortion might be 
expected when comparing them with vibration data. There are other 
data available in reports, which .are referenced in the published 
papers, but it does not appear that they would add much to what 


has been published. 



The nature of the published data is sunimarized in Table 1. 

From the table, we see that there are empirical formulae given in 
References 3, 5, and 9. The highest frequency is 5 k Hz, the 
largest SPL value is 150 db, and the highest temperature is 130® C. 
Thus, the data themselves are not directly useful for 15 k Hz 
soimd fields at 170 db and temperatures up to 700° C. Nevertheless 
we can compare our data with the empirical formulae, parametrize 
our data, extrapolate curves, etc., and consider what inferences 
might be appropriate for understanding the physical processes 
involved. 


The empirical formulae in Refs. 3, 5, and 9 involve the 
following dimensionless numbers: the Grashof Number Gr, the 

Prantl number Pr, the vibration Mach number Ma, the vibration 


Reynolds number the Nusselt number Nu. 


The definitions 


used in the papers are 


Gr = g D' 


BA t 


where B is the volume coefficient of thermal expansion (for an 
1 

ideal gas , B = -y ) 5 


Pr 


y Cp 
m 

k 

m 


where is the viscosity and Cp is the specific heat at 

constant pressure; 

Ma = _s 
C 

where is the velocity of acoustic streaming, arbitrarily 

chosen to be given by = ao, a being the maximum ' 
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TABLE 1 


Reference 

Method 

Frequency 
Range (Hz) 

’SPL 

Range (db) 

2 

Hor. St. Wave 

120 

110 - 117 

3 

Hor. St. Wave 

1100 - 4900 

0 - 151 

4 

Hor. Tr. Wave 

1000 - 500 

0 - 148 

5 

Vert. Vibrations 

54 - 225 

0 - 150 

6 

Hor. St. Wave 

1500 

0 - 146 

7 

Hor. St. Wave & Vib. 

100 

136 - 150 

8 

Hor. St.‘ Wave 

680 - 1090 

130 - 140 

9 

Hor. St. Wave 

645 - 672 

131 - 148 

10 

Vert. St. Wave 

710 - 1470 

125 - 140 

11 

Hor. & Vert. St. 
Waves 

710 - 1470 

125 - 140 

12 


0 

1 

20 - 40 

This Report 

Vertical St. 

15,000 

0 - 180 


Waves 


A/2D 

Range 

Temperature 
Range (°C) 

300 

49 - 52 

2-8 

20 - 130 

2-9 

100 

35 - 145 

35 - 125 

6 

115 - 165 

78 

50 

8-13 

30 - 50 

13 - 14 

85 - 93 

1.2 - 2.5 


1.2 - 2.5 



4.5 20 - 700 


Comments 

Temperature vs. angle 

Pictures ; graphs , 
formulae 

Pictures; Ah vs. SPL 

Graphs, formulae 

Local effects 

Compare vibration & 
sound effects 

CTl 

Local effects 
L vs. At; formulae 
Local effects 
Local effects 

Local effects 


OF POOR QUALITY - 


sotmd-particle amplitude of oscillation and uJ the angular 
frequency; 

V_D 

(Ref. 5) 

=-■ m 
V 'D 

Re = — ^ (Ref. 9) 

s V 

m 

with V^' = (2 TT a 0 ) <t>> ‘I’ being determined empirically at 

each temperature; and 
h 

Nu = AAt 
as before. 

There are two possible ways to define At. Since the sound 
field heats the coil even with no current supplied to it, the tem- 
perature t at Q = o will be higher than the ambient temperature 
o 

t . The difference is about 5°C for SPL = 170 and 12° C for 
a 

SPL = 180 db . The Grashof number would normally be defined so 

that Gr = 0 when Q = 0. On the other hand, the forced air flow, 

by acoustic streaming, involves air at ambient temperature (outside 

the standing-wave sound field) , so the Reynolds number would 

normally be defined for a mean temperature relative to ambient. 

The difference between t and t at 170 - 180 db is not trivial 

o a 

when one compares data plots using each one. If h is linear 
in one plot, it will not be linear in the other, even within 
experimental error. 

There is also an ambiguity in the definition of the Grashof 
number. Ford and Kaye (3,5) state that it is evaluated at the 
mean film temperature, presumably implying 8 = 
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uses 8 = 1/T . For At < 100“ C it made little difference, but 

a 

Gr peaks at A t ~ 250“ C using, Yuge's definition and at a lower 

temperature with a lower value using the definition of Ford and 
Kaye. 

The empirical formulae given in Ref. 3 are: 

h^ = 0.245 (At/D)^^"^ 

h^ = [h^ (At D^)™/D] E(a f)^ F]*^ 

1/4 

Nu^ = 0.485 (Gr Pr) 

Nu = B (Gr Pr)^^^ (Ma^F)^ 

V ^ 

Here the subscript 0 implies SPL = 0, V implies SPL =|= h and 
2 

h are in Btu/ft • hr • “F; b and B vary with frequency or, 

V VV 

more appropriately, with A/2D; m and n vary with A/2D; 
f = 2 ir co; and F is a geometrical factor included to compensate 
for the different values of the amplitude a over the surface of 
the cylinder, so it, too, varies with A/2D. The values of b^ 
and are tabulated and graphed, and the values of F are 

graphed. 

The formulae in Ref. 5 are: 


h = 
o 

0.255 (At/D)^^^ 

h 

b (At/D)°'^ (af) 

V 

V 

Nu = 
o 

1/4 

0.495 (Gr Pr) ^ 

II 

B (Gr Pr)®'^ Re 

V V 


n 


2 

with h and h also in Btu/ft • hr • °F. Thus the param- 
o V 

eterization used by Ford and Kaye for horizontal (3) and vertical 

(5) standing waves is similar, except for a change from Ma to 

Re . However, the number n is about 1/3, which means that 
5 

2/3 2/3 

there is a change from (af) to (af) and (Ma) to Re^ 

in going from horizontal to vertical. Note that the equations 

for h and Nu correspond to identical situations in the 
o o 

two cases. In both cases, the v-subscript equations are for 
"fully developed vortex flows" which, according to the authors, 
means SPL ^ 146 db. 

The parameterization of Ref. 9 is more complicated. Lee 
and Richardson (9) arrive at a formula: 

Nu/Gr^^'^ (1 + 1.61 = 0.372 [ 1 + A^(Re^')^/Gr]^'^'^ 

where Re (discussed above) is not calculated directly, but is 
s 

replaced by means of the relation 
= (0.372)^ A^0 ^ d^/v^ Gr 

with B^. determined independently at each temperature. Note 

2 

that the definition of Re gives Re (af) while Ref. 5 

s s 

has Re « (af ) . Since the dependence of Re on (af ) is never 
V s 

used, however, no comparison with Ref. 5 can be made. Lee 
and Richardson (9) claim to fit the data and Peebles (7) as 
well as their own, and demonstrate the fit with a graph. In 
Refs. 3 and 5, there is no demonstration of the degree to which 
the equations fit the data, except for h^. Since the data are 
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shown on a graph of log h is log At, it would be possible to 

investigate the quality of fit for (which I expect is rather 

poor) but I have not done so. 

The data of this report are tabulated in my laboratory 

notebook No. 2, and represented in Figure 1, which is a plot 

of log h vs. log At for the data of Ford and Kaye (3,5) is 

sho^m., and we can compare data in this range for SPL = 130, 140, 

and 150 db. (Note : the data at 130 db are so close to those 

at 140 db that they are not plotted) . The range of Grashof 

4 4 

numbers in Refs. 3 and 5 is 1 X 10 to 5 X 10 , while my 

range J.S 40 to 115; tjie range of Reynolds nubmers is 0 to 

3500 for Ref. 5, 0 to 300 for my data (up to 150 db) , but 

1/2 

typical ratios Re^/(Gr) are of the same order of magnitude. 

A distinct difference is that at higher values of At the 
Grashof number peaks, then decreases. 

An analysis of ordinary forced and mixed convection indicates 
(13) that one could expect to get 

Nu = f (Gr, Pr, Re^) 

with Gr and Pr appearing as a product. Consequently, I 

calculated Gr, Pr, and Re^ for my data. In the calculation, 

I chose Yuge’s definition of Gr ( 3 = l/l ) ^nd defined At 

relative to t . In Re , v = (acJ) is defined at t . By 
a v s a 

extrapolating to corresponding values of At at different SPL, 

I plotted log Nu vs. log Re^ at constant (Gr Pr) , as shown 
in Fig. 2. Although (Gr Pr)' is double valued, the curves below 
and above the peaks are clearly separated by different values 
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of Nu. Note that the data are very close to straight lines, 
and that the slopes are all close to the same value. For At 

< 150° Cj the slopes are about 0.56. For A t > 300° C, they 

are about 0.52. Fig. 3 shows a plot of (Gr Pr) is At and 

Re^ is At. The dependence on Re^ may be compared to the 
data for ordinary forced convection (See Ref. 13, p. 242). For 
Re from 0 to 400,000, the data are well fit by 

Nu = 0.43 + C (Re)™ with C=0.48, m=0.5 

for 0 < Re < 4,000 and c = 0.174, m = 0.618 for 4,000 < Re 

< 40,000. It is not possible to plot Nu vx. (Gr Pr) at 

constant Re- , because there are no common values of Re . 

V V 

0.52 

Consequently, I plotted (Nu/Re^ ’ ) vs. (Gr Pr) and got 

the graph shown in Fig. 4. The low temperature data would 

0 56 

obviously agree better if I used Re^ ’ ; nevertheless, the 

trend is clear. For At < 100° C we get 

~ 0.03 /Rg GrPr 

V 

while for At > 200° C we have 
Nu ~ 54 -/Re / (Gr Pr) . 

V 

Between 100° C and 200° C there is obviously some sort of 
transition. 

D . Further Work 

The figures included in this report contain no direct indication 
of agreement or disagreement with the published formulae. 
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However, since the low temperature dependence of Hu is 
clearly 

Nu <== (Gr Pr) 

there is no need to look for a comparison with the formulae of 
Ref. 5 where Nu «= Re . There is also little hope for 

V 

2/3 

Nu Ma as in Ref. 3, but this should be checked, as should 
the complicated formula in Ref. 9. It would also be worthwhile 
to investigate a parameterization like 

Nu = (Const.) +(Const.) (GrPr)™ (Re^)^ 

by anology with the good data on forced convection in Ref. 13. 

All of this can be accomplished fairly easily by computerizing 
the data and using a minimizing program which I have available. 

A good check would be accomplished by taking additional 
data at 145, 155, and 165 db. The data at 178 db are suspect 
because we don’t have a good measurement of SPL, but it would 
be worthwhile using the low-pass filter on the B & K amplifier 
and using the resulting SPL value as the correct one, just to 
see how the results would agree with data at 140, 150, and 160 
db. The 170 db data are also suspect for purposes of comparison 
because the field is not close to a sine wave and we are assum- 
ing a sine wave in our parameterization. 

Once the matters above have been Investigated or accom- 
plished, I plan to write a paper for publication (authors: Rush 
and Dean) in J. Acoust. Soc. Am. or J. Heat Transfer. I sug- 
gest the former, because they are more likely to be interested 
in heat transfer at high frequencies. I will, of course, try 
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to get a theoretical basis for the change in Nu vs. (Gr Pr) 
at 150® C, and for the differences between our results and others, 
before submitting the paper. 

Following the above work, which should take about one month 
or perhaps during the process of the above) I plan to investi- 
gate the stability of objects in the acoustic system at room 
temperature, introducing specific perturbations such as air 
currents and reflector changes, and quantifying the observations 
as much as possible. 
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Notation 


a = amplitude of sound vibration 

A = surface area 

A = cross-sectional area 
«) 

c = speed of sound (at t^) 

Cp = specific beat of air at constant pressure 

D = diameter of cylinder 

f = frequency (sound or vibration) 

g = acceleration due to gravity 

3 

Gr = Grashof ntimber = g D gA t/v 

m 

h = heat /-transfer coefficient = Q/AAt 

k = thermal conductivity 

L = length of segment of lead wire 

^ = vibration mach number = 

Nu = Nusselt number = hD/k 

m 

Pr = Prandtl number = u Cp/k 

m m 

P - root-mean-square sound pressure 
rms 

Q = convective heat transfer rate 

T 

Re^ = streaming Reynolds number == 

Re = vibration Reynolds nvimber = . V D/v 
V s m 

SPL = sound pressure level in db (re 2 X 10 dyne/ cm ) 

T = absolute temperature, K 
t = temperature, °C 

= velocity amplitude = au) 

2 

V^‘ = streaming speed = (af) 0 (Ref. 9) 

g = volume coefficient of thermal expansion = 1/T for ideal gas 
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At = temperatrue difference 


X = wavelength 
y = viscosity 

V = kinematic viscosity = u/p 

2 

0 = function of f, D, and V such that (af) 0 is a characteristic 

streaming velocity (Ref. 9) 
p = density of air 
0 ) = angular frequency = 2 tt f 

Subscripts 

0 implies Q = o (this report) or SPL = o (Refs. 3,5) 
m implies mean film temperature, average of surface temperature and 
ambient temperature 
a implies ambient temperature 

m refers to Kanthal wire leads to heating coil 
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APPENDIX B 


PROPERTIES OF A CONSTRICTED-TUBE AIR-FLOW LEVITATOR 
J. E. RUSH,* W. K. STEPHENS,* AND E. C. ETHRIDGE** 

*Department of Physics, University of Alabama in Huntsville, Huntsville, 
Alabama, USA 

**Space Sciences Lab, Marshall Space Flight Center, Alabama, USA 


ABSTRACT 

A constricted-tube gas flow levitator first developed 
by Berge, Oran, and Theiss shows promise both as a space- 
positioning device and as a levitator for ground-based 
work. We present results of laboratory studies which were 
designed to predict the behavior of the device in a low-g 
environment . 

0R1G5NAL IB 

INTRODUCTION OF POOR QUALITY 

There is much interest in levitation techniques which can be used in ground 
based research to study a number of phenomena. The effective levitation of 
liquid nonconductors remains an elusive goal in materials science. A 
constricted-tube gas flow levitator has been developed at the Marshall Space 
Flight Center by Berge, Oran, and Theiss [1]. Its advantages are that it is 
a simple levitator which is essentially orientation and gravity independent, 
will operate over a broad temperature range, and can be used to process both 
conducting and nonconducting materials. Solid spherical samples of a number 
of different densities at 1200°C have been successfully levitated. We have 
continued to study the properties of such levitators as a possible solution 
to this goal and as a possible levitator for low g and report here on our 
work to date (Sept., 1981), 

The levitator consists of a constricted (quartz) tube fed at one end by a 
source of heated air or gas. A spherical sample is positioned by the air 
stream on the downstream side of the constriction, where it can be melted 
and resolidified without touching the tube. The primary source of heat for 
the sample is the air itself, although secondary sources are also being 
investigated. The air is heated by being passed through a furnace, by being 
blown past a torch, or both.. 

The behavior of spheres in flowing fluid in a Ig environment has been 
studied by many investigators. The earliest work was on spheres in free 
flow [2], where the drag coefficient vs. Reynolds number Re was determined. 
Also investigated was the pressure distribution around the sphere and the 
behavior of the fluid as a function of Re. A more complex situation is flow 
past spheres in vertical, slightly tapered tubes forming the basis of flow- 
meter technology. Quite a bit of work has also been done on spheres in 
straight tubes and pipes, both experimentally [3] and theoretically [4]. 
However, the only published work on spheres in duffusers which we have found 
is by Schmidt and Springer [5j. 

In a straight, vertical tube one can maintain an unstable equilibrium in 
the vertical direction by balancing the gravitational force with a drag force. 
A relatively stable equilibrium exists laterally because of a Bernoulli effect 
due to the radial variation in fluid velocity, coupled w^ith a weaker Magnus 
force due to rotation of the sphere [6]. With a slightly tapered tube, as in 
a flowmeter, one c;,an easily produce stable equilibrium, but levitation is 
gravity-dependent. The effect of a constriction is to supply an upstream 
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(Bernoulli) force which can give stable equilibrium with or without gravity. 

The purpose of the work presented here is to study the properties of the 
levitator in order to predict its behavior in a gravity-free environment and 
at elevated temperatures. 

EXPERIMENTAL PROCEDURE 

We define the inner diameter of the tube as the constriction diameter 
as D 2 , and the diameter of the spherical sample as d. The equilibrium position 
of the sample in Ig and the stable pressure were measured at each end of the 
tube as a function of D 2 /D-j^, shape of the constriction, flowrate through 

the tube, and shape of the diffuser. We also studied the net axial force on 
samples suspended by a wire from an electrobalance as a function of position 
of sample and of the variables given above. Thus we were able to map the 
force in the neighborhood of the expected low-g equilibrium position. 

Since we are ultimately interested in levitating and melting samples at 
elevated temperatures we' want to know how the levitation forces are affected 
by temperature changes. For a given tube and sample size, with the system in 
mechanical and thermal equilibrium, the net .forces on a sample at any given 
point in the tube should depend only on the value of Re. By varying the fluid 
flow rate we can vary Re thereby modeling the effects of varying temperature 
on Re. Thus temperature effects can be studied by varying the flow rate. 

The tubes used in these studies had an inner diameter of approximately 
0.65 cm with ratios D 2 /Dj^ of 1/4 to 1/2. The ratios d/Dj ranged from 1/2 to 
7/8. For the shape of .the constriction, we found that a simple but adequate 
quantitative measure was the approximate angle a of the inside surface with 
respect to the tube axis near the equilibrium position of the sample. The 
values of a ranged from 7° to 18°. 

Each end of the constricted tube was attached to copper tubing by means of 
a plastic heat-shrinkable tubing. Fittings were mounted on each copper tube 
for flexible-tubing connections to pressure gauges and to a flowmeter on the 
upstream side. Gauge pressures were measured with and without samples inserted 
for the range of flow rates used (3 to 17 liters per minute) . At the lower 
flow rates we also measured differential pressures. We also measured the 
variation of pressure with sample position for suspended samples. The 
pressure at the downstream end was never significantly different from 
atmospheric pressure (AP< 0.01 psi) . The precise position of the sample was 
read with a cathetometer . 

The force measurements were made using a Cahn electrobalance. A fine wire 
was suspended from the electrobalance, passing through the tube, and attached 
to the spherical sample by means of a small bit of epoxy. In these measure- 
ment’s the sample material was not important, so we used assorted steel ball 
bearings . 

In each case the suspended samples were balanced before the airflow was 
begun, so that the forces measured were due entirely to the airflow. The 
suspended scimples, of course, could not rotate. Rotation would probably occur 
for most freely— positioned samples, but this effect would not greatly modify 
the restoring force (6) and would in fact improve stability. 

EXPERIMENTAL RESULTS 

In Table 1 we shov: some basic data for five constricted tubes which were 
used in these studies. 
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TABLE I 

Basic parameters of constricted tubes 
Inner 

Tube Diameter Constricted Ratio 

Number D^ (cm) Diameter D 2 (cm) D 2 /Dj^ Angle a 

1 0.650 0.256 ■■ O 9 IP 

2 0.650 0.350 0.54 9-5° 

3 0.650 0.167 0.26 7° 

4 0.650 0.157 0.24 18° 

5 0.676 0.157 0.23 17.5° 

A plot of the force, F, on a sphere In tube #5 is shown in Figure 1. Since 
the sphere is balanced, the force is due entirely to the fluid flow. This 
force is plotted against displacement of the sphere along the axis of the tube 
where the zero point of displacement corresponds to the point where the sphere 
touches the constriction and blocks the flow. The sphere used for Figure 1 
had a diameter of 0,475 cm (d/D2_=0.23) and the flow rate was 13.2 liters per 
minute, corresponding to an average velocity of 612 cm/s at diameter 
(upstream) . 



0.5 1.0 1.5 

Displacement (cj^) 


Fig. 1. Fluid force vs. displacement for balanced sphere. 
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Since ve want stable equilibritim, we are looking for a force vs. displace- 
ment curve with a large negative slope and a broad range of positive/negative 
values to the left/right of zero. We also want the maximum and minimum values 
of the force to be reasonably large in comparison to typical sample weights. 

For a sample of the size used in Figure 1, which has a good d/D^ ratio for 
stability, the weight for the steel sample is 440 mg. Thus the negative value 
for the F (10^ mg) is quite acceptable. Since the sphere is easily levitated 
at this flow rate, the positive value is acceptable also. The least desirable 
feature of Figure 1 is the small value of the equilibriiim displacement, but 
this can be adjusted by varying the shape of the tube. 

-In Table II we summarize* the results with the various tubes in terms of 
varying diameters of samples. In each case the flow rate chosen for comparison 
was based on the minimum value of the force for the range of flow rates used. 


TABLE II 

Best results for force minimum with each tube and sample 


Tube 

Number 

3-) 

Sample 

(Sphere) 

Force 

Minimum (mg) 

Flow 

Rate (2/m) 

Equilibrium 
Position (mm) 

1 

A 

None 




B 

-380 

14 



C 

None 



2 

A 

None 




B 

-70 

5 



C 

None 



3 

A 

-40 

4 

2.0 


B 

-620 

9 

2.5 


C 

-450 

14 

2.5 

4 

A 

-820 

14 

1.3 


B 

-960 

13 

0.8 


C 

-990 

13 

0.5 

5 

A 

-990 

14 

0.17 


B 

-1020 

13 

0.7. 


C 

-1020 

13 

0.6 

Sample 

diameters: A, 

0.600 cm; B, 0.475 cm; 

C, 0.356 cm 



One is also interested in the motion of the sample in the tube, primarily 
to avoid contamination by touching the walls and to avoid liquid sample break 
up. For the better diameter ratios obtained from Table II, there was no 
apparent motion of the sphere at equilibrium. 

In Figure 2 we give some typical results for a freely levitated solid sphere 
at equilibriimi. The larger displacements occur, of course, when the tube is 
inverted. From the small displacement on inversion (see Fig. 1) it is clear 
that One can choose flow rates such that the gravitational force becomes in- 
significant, as was demonstrated by Berge, Oran, and Theiss. 

CONCLUSIONS 

From the results given above, we conclude that the constricted— tube levi— 
tator can be used successfully as a positioning device for solid spherical 
samples at low g. 

For operation in thermal equilibrium at high temperatures, we note that the 
important fluid parameter is the kinematic viscosity. If air is heated from 
room temperature to, say i200“C, the kinematic viscosity increases by a factor 
of 14, To maintain a given value of the Reynolds number, the flow rate would 
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have to be increased by the same factor, for a specific geometry of tube and 
sample. Thus, to maintain stable equilibrium, one increases the flow rate as 
the air or other gas is heated. The feasibility of this process has already 
been demonstrated in 1 g by Berge, Oran, and Theiss. 



4.0 $D 1X0 16.0 

Flow Rate. (Ljr£a5//vuN) 


Fig. 2. Position of freely levitated steel sphere vs. flow rate at 1 g. 

The other stability problem which must be considered for processing of 
samples is the change in shape of the spherical sample as It melts. The solu- 
tion to this problem involves selecting a shape for the constriction so that 
the solid sample does not spin too rapidly and does not contact the tube on 
meltinp- From the data of Tables I and II, one can see that by decreasing 
the C 07 Etriction angle (changing from Tube 5 to Tube 3) one can maintain a 
sufficiently large force and significantly increase the separation between 
tube and sample. We are currently beginning a study of the stability of 
melted samples at 1 g. 
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ORIGINAL RAGE gS 
OF POOR QUALiTY 
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RAD I AT, RADI AT 
08/31/'82-09:1 2 (, 0) 

...4-^ REA L * S . D.E , H M D E I 

2. KEAL*i, DTA 

tX.TEH4.AL . f!^’DE 1 

CO M n U K 1 M / _ X H E; T £, ( 

COMMON /I NT/ 


3. 

4. 

5 . 

6 > 

8 . 


23, 

24. 


26, 
27, 
..Z8 

29. 

30. 
34- 
3 2-^ 

.3i»- 


-3-4.^ 
35 . 
36-. 
-33.- 


1 

..J. 

1 


' » 

PSI ( 50, 50) 
, ST(50) 


■ 1 ) 


, THETAO.C 50) 
, PSI0(50) 


lAU ( 5,0J.- 


,.x:( 5 0.) 

, El noo) 


-I (50) 


,TSTAR(50) 
, U-( 5 0-, 5 0 ) . 


TPRI mE ( 5 0) 


9 . 

1 0 , 

read ( 5 , 101 ) NX, NT, M 

liLl EJ-tAJi: ( 

1 1 .■ 


READ ( 5 , 102 ) AN, TAUO, 

12 . 

102 

FORMAT ( 5 FIO.O) 

-4 3 ...— 


-w-S-tJ-E ( d.,J. 5 -. 1 J.. A.T, -.TAiiCL,- 

1 4 . 

151 

FORMAT ( / / '♦HN = , F 5 . 2 ,4 

1 5 . 


OELX = 1 ./NX 

—1 O ' — • 


p F I , T A Lf =._.T A U 0 * a ELX . 

17 . 


R = DELtAT/ D bLX *-*2 

12 . 


NXP = NX + 1 

.- 19 -._- 

■ - - 

. kX.PP= NX +.2 

20 , 


vO 201 1 = 1 , NXPP 

21 . 


X(I) = (I- 1 )*DELX 


DELTAT, ALPHA, 9ETa 


DEL-TAX. 


F5 . 2 , 4X , 1 OHOELTA T = ,E8.2//) 


-)- 


193 


— XALLC i ). A U Q-tX C . 

WRITE (6,193) TiU(I) 
FORMAT (60X,E8.2) 

-GAMtXlA .=.-.0.31... 

V.= - GAMKA^X(i) 

U Cl , 1 ) = EXP (V) *X ( I ) 

XmLIMUE 

NTP = NT + 1 
00 211 J=1,NTP 

._I_( J.)„ = -J.J ^1 ) * D E L.T A T 

-TSTA.RCJ) = ALPHA’^TCj.) 

_T RRJ M E-C j. ) - .=.. b.LT A.* T 5 .lAiU J ). 


-2_U_CDiLLLNiJX. 


NXD = 2.* NX + 4 
DO 221 ]=2,NX0 
_XJ._^. .(_L-J.) * D £.LX. . 


33. 


TA = TAUO*XA 

39. - 


DTA = TA 

40 . 


D . . = AIM D Ei .(2 , OJ.A , 1 E.R.) . ... 

41 . 


.IF (IER.ME.130) GO TO 310 

42. 


WRITE (0,131) 

43 

.18.1 

.FORMAT <//.//5H lERO) 

4 4. 

310 

continue 

45. 


IF ( I ER. NE. 1 31 ) GO TO 311 

46-. 


-U'RJLT E . ( 6 , .1.ii.2 ) _ 

47. 

182 

FORMAT (/ ///CH09ERF LOW) 

48. 

31 1 

CONTINUE 

49. .... 

. ^ 

El(I) = DE 

50. 


WRITE (6,191) EKI) 

31 . 

1 91 

FORMAT (30X,E8.2) 

52. 

22 . 1 . 

_Ci).NJ.iriUE. .. - _ 

53. 


DO 241 J=1 , NTP 

> 4 , 


DO 251 1=2, NXP 

5-5.. 


-..T.HEXA (1,J0 =.U(L,J)/X(1) .. 
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56. 251 COMIWUE 

57. FN' = 0.0 

. S.a,. HO 261 1 = 2, NKP. 

59. DO 271 K = 1 , NXP 

60. 1H6 = THETA (-L,J ) •* 

o1 kA„= lAFiSd-k) . 

62. kB = I + K 

65. IF (KA.Fifc.O) GO TO 253 

6A-. £I-C!CA) = Q-. Q. 

o5 . - 2 53- CONTINUE 

66. EDIF = EI(kA) - LJ(kB) 

_6-7-. f-C-N- = - 2 . * T.HiL>' ED i-F 

68. ■ - - IF (K.EQ. 1) GO TO 252 

69. IF (K.EQ. NX) GO 10 252 

...7.0- F-CN— =_3-. * F.CM ...... . . .... 

71, . 252-CONTlNUE 

7 2. F N = F r-f + F C N 

7 3.. 6-Z-1 — Caid X-(\ U-L- . - - . - _ 

7A. TG = 0 . 375*0ELTAU-* F N 

75 . , 0(1 , J+1 ) = 0. 0 

76. 6 U.(NXEP , j-t 1J. = U(NX,jn)-- . 

77 . TH4 = 1 hETA ( I , J ) * ^4 

78. PHI = -4. *TAO ( I ) * TH4 + TG 

. -_Z9. ^LLU = _XAU 0. /.At'Or-P.hJ 

80. WRITE (6,192) FLUX 

SI. . 192 FORMAT (50X,ES.2) 

. _a2-. UE_=L, UjL l±.\ , JJ . , . . _J 

S3. ur. = U ( 1 - 1 , J } 

ii4 . UN = U ( I , J ) 

_ Si. 0-(.L., J+JJ = .r6*_LU.P--UEl). .( 1 . -2 ..*R,0 * UN..+ . 0 EL.T A„X * F LLLX. 

86. 261 CONTINUE 

o7. U(1,J + 1) = ^.^OCEiJ-tl) 3.-*U(3,J+1) + U(4,J + 1) 

—Si. 2-4J— CQNJ.iNUE. .. . 

89. WRITE (6,.1oCO 

-90. .t6(l-F..0kMAT (/.///,) 

2J-. W.RITE (j6,T6.1)- 

92. '161 FORMAT ( 2>' , 1 HX , 1 OX , 1 HT , 1 1 X , 1 HU , SX , 5HTHETA/ ) 

93. '' - DO 301 J = 1 , NT 

2_4- D-Q_3 0i..l=-1 , MX . 

95. --- WRITE (6,1o2) X(I), T(J), U(I,J), THETA(1,J) 

96. .. 162 FOKHAT ( F 5 . 2 , 4X , E8 . 2 , 4X , E8 . 2 , 4X , El 2 . 6) 

. -97- 302. CON-TiNUE .. - . , . 

98. 301 CONTINUE 

99. STOP 

1 C 0 END,. . 



